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ABSTRACT 

The accuracy with which a troblle lunar laser station can be 
located is the subject of a modelling study presented here. The 
influence of the numljer and accuracy of fixed lunar ranging stations, 
the uncertainty in polar motion, and data loss due to weather and 
similar factors has been considered, and the results are given in 
cartographic form. In general, all three coordinates (for coordinates 
to latitude ± 60 ") are determined to better than the pole uncertainty, given 
three nr more fixed sites and reasonable .eather. This result indicates 
that one or more mobile stations will be extremely suitable for the 
study of geotectonics. 
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1. INTRODUCTION 

Until recently, the lunar laser rankling project has concentrated 
on the improvement of the astronomical model: the lunar orbit, lunar 

rotation, and coordinates of the reflectors and of McDonald Observatory. 
Refinements in the physical model, particularly the librations, now 
permit the fitting of 5 years' observations to a contimiocJ orbit with 
a mean residual of about 3 nsec (45 cm one-way range), the observations 
having formal uncertainties between 0.6 and 2.8 nsec 111. 'ITiese 
refinements have permitted also the first attempts to determine the 
variation in UTO at McDonald [2]. 

With the first lunar-laser determination of an intercontinental 
baseline [3], the commencement of regular operation of a fixed station 
in Hawaii, the advancement of fixed stations in Australia and France, 
and the soon- expected authorization of a mobile station, the lunar 
laser project has entered a new phase. This technique is on the verge 
of becoming a tool foi global geophysics. In support of this aspect 
of the program, we have conducted a limited sensitivity study to 
estimate the formal precision with which a mobile station can find 
its own geocentric coordinates under a reasonable observing schedule. 

The problem of modelling has been conceived in two parts, whlcl. 
we designate as geometric and statistical. The geometric is, quite 
simply, the uncertainty due to Incomplete knowledge of the orientation 
of the terrestrial coordinate system and is essentially unrelated to 
the conditions at the mobile site; for the present study, we have 
included only the polar motion uncertainties in the geometric problem. 
The statistical uncertaln+'y is provided by a covariance analysis of the 



system of normal equations for the solution of the mobile station 
coordinates from an hypothetical observation sot. Of course, the 
validity of such a process is at least in part dominated by the 
assunptlons that go into it. We have tried to adopt a r(*asonable 
observing schedule and to account for the vagaries of weather in a way 
consistent with experience at McDonald, but more can be done in this 
direction. 

2. THE GEOMETRIC PROBLEM 

In the present study, we have considered two levels of knowledge 
for the polar motion, one in which the uncertainty is ±50 cm, which is 
somewhat larger than the currently probable BIH errors, and the 
second with ±5 cm uncertainty, a pessimistic estimate of the level 
to be expected from pole determinations based on a network of lunar 
ranging stations. In addition, it has been supposed that each fixed 
station is capable of determining its astronomical variation in latitude 
to the level of Internal systematic errors (± 1 cm/sinrp for a 3- cm 
station). The justification for this supposition is in the averaging 
to be- cxptTienced over each pass of 6-10 hours, combined with the 
likelihood of correlation times in excess of one day for this 
phenomenon. The effect of this assumption is that the latitude of any 
point is better determinable along the meridian of a fixed station 
than would be permissible elsewhere. To see this, we will adopt temporarily 
an Earth- fixed coordinate system based on the true equator and the 
meridian of tlit' fixed station. The BIH equations for th(> effects of 
pole displacements x and y are 
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t\ = - {X sin X + y cos X } tan re 
e e p ' 

( 1 ) 

8<p « + {X cos Xp - y sin X^} 


where X^ and ?p are east longitude and latitude, respectively. A3thmig)i 
we suppose that the uncertainties 5x * 6y = Y, tl.'> aJditional infm>nKi- 
tion in the meridian of the fixed station makes it desirable to 
write the BIH-like system for the mobile site: 


?»X„ = - {x"' sin A + y’ cos A } tan m 
m m m m 


( 2 ) 


where A * X -X.., subscripts m and f denote mobile and fixed sites, 
x' and y' are the components of pole position in the new frame. Thus, 
if X - 1 1/sin CP ! , the resulting? uncertainties are 


Uf eXm) = { |X sin A I + I Y cos i tan c::^ 
Uf (cPm) = { lx oos A I + I Y sin A I ) 


(3) 


The total uncertainty assigned to the geometric aspect of the problem 
is then computed as the harmonic mean of the individual fixed station 
contributions 
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If the effects of all fixed stations were Identical, this would 
coi’respond to reducing the uncertainties by the square root of the number 
of fixed stations. 

3. THE STATISTICAL PROBLEM 

To sufficient accuracy for this study, one can represent the 
topocentrlc distance p of a reflector on the lunar surface by 

p = R - o cos 6 cos H - Z sin 6 (5) 


where R is the projection on the line of sight of the geocentric 
reflector distance, a the telescope spin axis distance, z its equatorial 
distance, H the local hour angle of the reflector and 6 its topocentrlc 
declination. This relation readily gives the partial derivatives of 
topocentrlc distance with respect to telescope coordinates: 

•^ = ^ = a cos 6 sin H 


- cos 6 cos H 
bo 


M = _ 

bz 


sin 6 


( 6 ) 


Each observation gives an equation of condition relating the (observed 
minus computed) range residuals Ap to the required improvements , 

Ao, Az to the estimated telescope coordinates; 
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Hach observation contributes to the system of nomal equations In a 
classical Gaussian least squares analysis. 

a) Scheduling 

One of the features of the lunar motion that Is at once an 
advantage and a nuisance Is the fact of Its slowness In Its orbital 
motion. The advantages are related to the length of an observing 
passage, permitting the buildup of good statistics on the diurnal 
rotation of the Earth, which reflects Into good stability for solutions 
of station longitude and spin axis distance. At the same time, howc*ver 
it requires a long time to cover the full geometric range normal to 
the equator. For this reason. It is supposed that a mobile site will 
be occupied for a minimum of 24 days. Despite the fact that the mobile 
station will be a dedicated facility, we have conservatively adopted 
here the same daily schedule followed at McDonald: three ranging 

periods of 1 hour duration centered about liour angles -3h, Oh, +3h. 

This gives a total of 72 possible observations, but there are other 
boundary conditions that serve to reduce the number actually considered 

As a small step towards reality, the time of the "observation" 
is not fixed rigidly, but is situated within tlie available hour-’.'ide 
interval by a random number generator function. 

b) Zenith distance restrictic.i 

As a first restriction on the availability of an hypothetical 
observation, we adopt the geometric consti-aint imposed at McDonald by 
the Federal Aviation Administration: the laser may not be fired at 
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zenith distances greater than 70*. The official restriction may be 
less In many sites, but below this limit the air r^iss becomes so 
great that refraction modelling is a problem. The effect of this 
restriction is to reduce the available hour angle spread when t)ie 
Moon is in the opposite hemisphere from the station. Again trying to 
be conservative, we simply exclude the corresponding "observation," 
rather than rescheduling closer to the meridian. An example of the 
results of this restriction Is given in Appendix A. 

c) nffects of weather 

Optical techniques such as laser ranging are stopped by clouds 
and by atmospheric turbulence (the McDonald laser Is not fired If the 
seeing is worse than 6") . Our program accepts as Input a weather 
factor WFAC, which characterizes the incidence of usable sky. In 
using this factor, we try to in^^lude effects to account for both the 
medium-term randomness of weather with the fact of day-to-day correla- 
tion due to the size of weather patterns. This is accomplished by 
generating, for each observation window, a random number which is 
combined with a correlation function and compared with the weather 
factor. If the computed number is larger, the "observation" is 
nullified, contributing nothing to the system of normal equations, and 
weakening the statistical determination. The correlation function has 
two forms, chosen by experimentation and depending on the value of the 
weather factor. The derivation of this function is detailed in 
Appendix B, along with a series of large-scale sky-clearness maps tiiat 
are a byproduct of this study. 
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d) Determination of statistical uncertainty 

rinally, the effects of individual random numbers arc- reduced 
by averaging the res'ilts of ten computations. Admittedly, tills is not 
an Ideal means, but It seems to be a usable compromise between the 
desirable and the possible. 

Each of the available observation times that Is not eliminated 
by the zenith distance test or by the weather model contributes Its 
share to the noitnal system. The result Is a set of covariances that 
represent the determlnablllty (In the language of statistics, tlie Inverse 
weight ! of each of the coordinates, or their expected uncertainty, due 
to the effects modelled. 

4. RESULTS 

The final uncertainty. In this study. Is determined as the root 
mean square of the statistical and the geometric contributions. We 
present here a selection of cases that have been computed to Illustrate 
the principal res’ilts. In all cases, we assume 3-cm (0,2 nsec) 
capability for fixed and mobile stations. This is a factor of 3 better 
than current McDonald operation, but equals the Hawaii design specifica- 
tion; eventual operation is expected at this level at McDonald also. 

The variation between cases depends uniquely on the polar motion model 
and the number and location of fixed stations. The results are give. 
i.n graphical form, grouped by coordinate (radial distance, longitude, 
latitude, spin axis distance, equatorial distance). Noti* tliat i*adial 
distance Is a function only of weather factors, thus only two charts 
are given. For longitude and latitude, we present the follm^lng cases: 
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a) Worst case - 

This case assumes 50-cm polar motion information. All other effects 
are trivial in comparison. 

b) Statistical case, WFAC « 0. 5 - 

Assumes no polar motion uncertainty, which is unrealistic, but 
permits one to see the intrinsic level of the statistical contribution 
when of the observation periods are weathered out. 

c) Statistical case, WFAC ■ 1.0 - 

An unattainable best case, demonstrating the effect of a good weather 
site, relative to the preceding. 

d) Three fixe 1 itatlons (Texas, Hav.’all, Australia), WFAC *=0 5- 
The most probable configuration In the immediate future, with a 
believable weather factor. 

e) Three fixed stations, WFAC = 1.0 - 

f) Four fixed stations, IVFAC = 1.0 - 

Showing the gain to be realized from an additional station, in southern 
France. 

The last three cases are also illustrated in axial distance 
and equatorial distance. 

5. DISCUSSION 

The present study Is neither definitive nor without its flaws. 
In general, we have neglected all systematic errors related to the 
orbit and rotation of the Moon. A justification for this is tlie 
existence of tlie network of fixed statior.s, wJiich may be I’egarded as 
continuously monitoring these effects. The purpose of the mobile 
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station is puroly terrestrial, and one can consider that the present 
discussions concern only the coordinates of the mobile station relative 
to the systimi defined by the fixed network. We do not consider this 
to be a serious problems. Systematic errors related to mobile station 
hardware is more serious: wlien the calculated coordinate uncertainties 

arc* smaller than those expected from the systematic error biidget 
for the mobile station, then of course they should not be taken too 
seriously. Preliminary studies indicate that the estimate of 5 cm 
polar motion from a network of 3 cm fixed stations is reasonable, 
but more study is required in this area, and is indeed shortly to be 
undertaken. The use of a 50% weather factor is arbitrary, although 
it I’orresponds closely to the observation loss due to comparably 
random factors during several years’ operations at McDonald. 

Among the demonstrated advantages of lunar laser ranging for 
the study of secular variations in the orientation and distortion of 
the Earth is the combination of feasibility of dally operation over 
long periods of time and the extreme stability of the orbit of the Moon. 
It is not convincing to study long term secular changes by means of 
discontinuous models for the physical behavior of the observed 
object. It is precisely this point that makes our natural satellite 
a more suitable target for some purposes than any oth^r Earth satellite. 
Present data analyses using a single uncorrected, unrectifled, 
continuous orbit for the Moon over an Interval of more than 5 years 
give mean residuals less than 40 cm without the determination of ITl 
and polar motion from these data. The present study was undertaken 
within the context of that situation as a first step in evaluation 
of the potential utility of a mobile station for geodynamic studies. 
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The* lnte>rpretatlon of tht» prosent rosults may bo roRarded, as with 
all ('cvarlance studies, as primarily qualitative rather than quantitative. 
Nonetheless it seems clear from the results that a mot •' t lunar laser 
ranRinR station will be capablle of locating itself in all i c .ordinates 
to an accuracy comparable with its own observational uncertainty. The 
degree to which this can be achieved will depend, of course, on several 
factors not all of which concern the mobile site itself. In particular, 
it is absolutely essential that a high precision source of polar 
motion information be operational; we envision that this of course 
will be a network of fixed lunar laser ranging stations, some of which 
should be at relatively high latitude. The availability of good 
weather at the mobile site is Important, quite naturally, but the 
requirements are not extreme. 
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Appendix A: Effects of Zenith Distance Restriction 

For the purposes of this study, the motion of the Moon in 
declination was modelled by a simple sine curve. For each of the 72 
observing windows, the computed declination was combined with the 
lunar hour angle and the latitude of the mobile site to determine the 
zenith distance of the Moon at that instant. It Is obviously 
necessary that there be some geometric limit that will eliminate the 
possibility of making some observations. If there were no stricter 
limitations, this would be provided by the hu.-lzon; l.e. the Moon 
cannot be observed when its rienlth distance is greater than 90®. In 
fact, there are stricter limitations. They may be physical, such as 
the difficulty of making observations through the extended air mass 
near the horizon. They may be environmental, such as the 70® zenith 
distance imposed on laser operations at McDonald Observatory due to 
conslderatl' ns of aircraft safety. Taking this latter as our criterion, 
we have automatically rejected any potential observation less than 
20® above the horizon. This effect is obviously a function of the 
latitude of the mobile station, and it takes two general forms. 

Quite obviously, when the algebraic difference of the station latitude 
and lunar declination exceeds 70®, there is no possibility of an 
observation at any hour angle; that portion of the Moon's orbit is 
simply unavailable from that site. This will only occur for stations 
at latitudes higher than about ±40®. Even ;inder more favorable 
conditions, however, and over a much wider range of latitude, obser- 
vations may be geometrically possible at meridian passage but not at 
the +3 hour windows. Thus, while virtually no observations are lost 
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at the equator, rougly 30% are eliminated by this criterion at 
latitude 45“ , and fully half at latitude 60". A station at one of 
the poles operating under this geometric constraint, would obtain 
observations only when the Moon was at its extreme northern declination 
and on the meridian. Table I illustrates the situation for the two 
intermediate latitudes. Note that some effect of the random variation 
within the one-hour windows is noticeable in the entries for 6 = 45" 


and days 13, 23 and 24. 
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Table I. Effect of 70“ zenith distance limit at two latitvides 


H = -3^ H = 0 M = +3^' 



Appendix I'; Weather Modelling 


As Indicated in the main text, the actual number of observations 
entered into the formation of the system of normal equations is 
reduced as a function of a weather factor WFAC assigned to that 
particular site. The mechanism by which this is accomplished is to 
use the assigned weather factor together with a correlation function 
intended to simulate the day-to-day nonrandomness of major wf*ather 
patterns to obtain a specific weather index WT for that particular 
observation. W is computed as the difference between a function of 
WFAC and the normalized product of a random number generator. If 
its value is negative then the observation is rejected ns having been 
clouded out. Since WFAC has, by definition, maximum value unity, one 
would expect that the desired result of selecting (72 x V\TACj observa- 
tions could be achieved simply by defining WF as the difference be- 
tween WFAC and the random number at each step, 'lliis, however, would 
achieve no correlation from one observation window to the next, wliich 
is clearly an unrealistic situation. We have attempted to introiluce 
reasonably believable correlations on two levels, internal to a given 
day and from one day to the next. This required that the first obs(>r- 
vatlon window each day be treated differently from the succeeding ones. 
For these latter, it was found to be a very effective representation 
of the dally correlation to set 

WF = WFAC - X + W 

where X represents the random number and W the value of \\T from the 
previous observation window. The situation for the initial window 
of each day was considerably more difficult. Tlie simple scheme used 
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within a day, If extended to all observation windows, resulted In an 
unacceptatly high level of snowballing and thus very poor overall 
representation. After an extensive heuristic study. It was concluded 
that no one single technique would suffice for the entire range of 
weather factors. To produce reasonable results. It was necessary' to 
adopt two different functional relationships for WF to cover the 
range of interest for WFAC. For weather factors In the range 0.2 - 

0. 5, the following relation was used 

WF » WFAC - X + I + (0.25 - WFAC^) 
while for higher weather factors, this relation was adopted 

WF = WFAC - X + ^ + (0.23 - 

As bef ire, X represents the random number and W represents the value 
of WF from the previous observation window, which in this case was 
the previous day. None of the schemes tried were capable of representing 
the case WFAC = 0.1, but this Is not regarded as a serious drawback. 

Areas where the weather is so bad as that will be studiously avoided 
by a motile station. 

Sensitivity computations were performed over a wide range of 
values of the weather factor, and it was originally intended to produce 
maps based on large scale climatic records. It was pointed out, howe' 
that it was virtually always possible to find specific sites in any 
region where the local weather was considerably better than tlie 
regional average. As a result of discussions with the requestors of 
this study, it was decided that only the two unifoimi cases of WFAC = 

1. n and 0.5 be included the main body of this report. Nonetheless, 
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during the early stages of this study a certain amount of effort was 
expended on the compilation of global weather statistics. This 
Information has been reduced to maps of clear sky fraction (given - In 
tenths) over a 15* grid. These data have been compiled on a 
seasonal basis, and there Is Interesting information to be foiiml her«*. 
Consequently, even though they were not used In tin* preparation of 
the results given here, we close this appendix by displaying the 
four seasonal clear sky maps plus a fifth showing the annual means. 
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